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ABSTRACT 

Recent observations have constrained the galaxy ultra-violet (UV) luminosity function 
up to z ~ 10. However, these observations alone allow for a wide range of reionization 
scenarios due to uncertainties in the abundance of faint galaxies and the escape fraction 
of ionizing photons. We show that requiring continuity with post-reionization (z < 6) 
measurements, where the Lya forest provides a complete probe of the cosmological 
emissivity of ionizing photons, significantly reduces the permitted parameter space. 
Models that are simultaneously consistent with the measured UV luminosity function, 
the Thomson optical depth to the microwave background, and the Lya forest data 
require either: 1) extrapolation of the galaxy luminosity function down to very faint UV 
magnitudes ~ —10, corresponding roughly to the UV background suppression 

scale; 2) an increase of the escape fraction by a factor > 10 from z = 4 (where the best 
fit is 4%) to z = 9; or 3) more likely, a hybrid solution in which undetected galaxies 
contribute significantly and the escape fraction increases more modestly. Models in 
which star formation is strongly suppressed in low-mass, reionization-epoch haloes of 
mass up to Mh ~ 10^° M0 (e.g., owing to a metallicity dependence) are only allowed 
for extreme assumptions for the redshift evolution of the escape fraction. However, 
variants of such models in which the suppression mass is reduced (e.g., assuming an 
earlier or higher metallicity floor) are in better agreement with the data. Interestingly, 
concordance scenarios satisfying the available data predict a consistent redshift of 50% 
ionized fraction Zi.oion(50%) ^ 10. On the other hand, the duration of reionization is 
sensitive to the relative contribution of bright versus faint galaxies, with scenarios 
dominated by faint galaxies predicting a more extended reionization event. Scenarios 
relying too heavily on high-redshift dwarfs are disfavored by kinetic Sunyaev-Zeldovich 
measurements, which prefer a short reionization history. 

Key words: cosmology: theory ~ intergalactic medium - reionization - galaxies: 
high-redshift - galaxies: formation - galaxies: dwarfs 



1 INTRODUCTION 

The installation of the Wide Field Camera 3 (WFC-3) on 
the Hubble Space Telescope (HST) has recently improved the 
efficiency of searches for faint z > 7 galaxies by more than 



an order of magnitude (e.g., Bouwens et al.pOlOa McLure 
|et al.|[20T0} [Bunker et aL]|2010| ). As a result, deep WFC-3 
observations have provided new measurements of the rest 
frame ultra-violet (UV, ~1,500 A) galaxy luminosity func- 
tion at these redshifts. These measurements are particularly 
important since galaxies are the most likely sources of hy- 
drogen reionization (e.g., Madau et al. 1999; Faucher-Gigere 
et al. 2008a, b). Nevertheless, it is difficult to robustly trans- 
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late these measurements into predictions of the reionization 
history, because of significant uncertainties in the spectral 
energy distribution (SED) of the galaxies, the fraction of 
ionizing photons that escape into the intergalactic medium 
(IGM), and in the contribution of fainter, as of yet unde- 
tected galaxies. 

Because of these uncertainties, it has been unclear 
whether star-forming galaxies can actually reionize the Uni- 
verse by z ~ 6 (as required by the transmission of the 
Lya forest at lower redshifts; Fan et al. 2002; Becker et al. 
2007, although see McGreer et al. 2011) and account for the 
Thomson scattering optical depth to the microwave back- 
ground implied by the latest, 7-year Wilkinson Microwave 
Amsotropy Probe analysis (WMAP-7; |Komatsu et al.|2011[ ), 
corresponding to a redshift of instantaneous reionization 
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-Zroion = 10.4 ± 1.2[^Even if galaxies are in fact the dom- 
inant re-ionizing sources, it is not clear to what extent faint 
sources below the detection hmit of existing observations are 
needed. 

The amount of star formation taking place in low-mass 
dark matter haloes is not only relevant for reionization, but 
also for our understanding of galaxy formation and evolu- 
tion in general. Indeed, several lines of evidence suggest that 
star formation in such haloes is suppressed, at least in cer- 
tain regimes. For instance, it is well known at lower red- 
shifts that the baryonic mass fraction in low-mass haloes 
is strongly suppressed relative to f2b/f2m (e.g., Conroy & 



Wechsler 2009 Guo et al. 20101. This baryon deficiency 



is commonly attributed to a combination of feedback pro- 
cesses, such as galactic winds, and suppression by the photo- 
ionizing background (e.g ., Dckcl fc Silk 1986, Efstathiou, 
1992} [Murray et al.|200"5'[ [Faucher-Giguere et al.|2011[ ). Ob- 



servationally, there also appear to be far fewer dwarf galaxies 
in the haloes of the Milky Way and M31 than the number 
of dark matter sub-haloes capable of hosting them predicted 
in A^— body simulations (e.g., [Bullock et al.||2000| [Madau 



et al.|2008 l, suggesting that some process inhibited star for- 



mation in the dark sub-haloes. Recently, theoretical models 
have also suggested that star formation may be specifically 
suppressed in low-mass haloes at early times due to a metal- 



licity dependence of the star formation efficiency ( Robertson 



& Kravtsov 



Dekel[[2011 



[20081 [Gnedin fc Kravtsov[[2010[ [Krumholz fc 
Kuhlen et al.[[20T2 l. If star formation is indeed 



strongly suppressed in early dwarf galaxies, then it may not 
be possible to rely on them to reionize the Universe. It is 
thus necessary to clarify the importance of those galaxies 
for reionization. 

The primary goal of this paper is to examine the exist- 
ing observational constraints on hydrogen reionization and 
its sources, and to systematically determine which scenarios 
are (and are not) allowed by the data. A main distinction 
of our study relative to recent analyses (e.g., Bouwens et al. 
[2011c| [Shull et al.[[20TT| [Jaacks et aL][20T2l ) is the inclusion 



of lower-redshift Lya forest data (see also iMiralda-Escude 
2003; Bolton fc Hachnch 20 071 |Fauche r-Gigu&r e"et al.|2008a[ ' 
Pritchard et al.„2010, ,Haardt fc Madau„2011j ). The mean 
transmission of the Lya forest, which is set by a balance be- 
tween the ionizing background and recombinations, has the 
advantage of being a complete probe of the ionizing sources. 
The total instantaneous rate of injection of ionizing pho- 
tons into the IGM can be measured from the Lya forest 
without recourse to assumptions on the escape fraction or 
extrapolating the contribution of faint sources, two of the 
principal uncertainties affecting traditional analyses based 
on the galaxy UV luminosity function. We also include re- 
cent constraints on the duration of reionization from mea- 
surements of the kinetic Sunyaev-Zeldovich (kSZ) effect by 
the South Pole Telescop^{SPT) high-resolution microwave 
background experiment (Zahn et al. 2011; for a recent pa- 
rameter space study of the kSZ signal from patchy reioniza- 
tion, see Mesinger et al. 2011). 



^ In reality, the epoch of reionization is expected to be extended 
in time and the Thomson scattering optical depth only provides 
an integral constraint on reionization. 
^ http://pole.uchicago.edu 



While the Lya forest data are mostly restricted to 
z ^ 6, when reionization is probably complete]^ they pro- 
vide valuable constraints in two ways. First, realistic reion- 
ization scenarios should continuously connect to the post- 
reionization IGM probed by the forest. Second, measure- 
ments of the galaxy UV luminosity function (analogous to 
those directly probing the epoch of reionization) are avail- 
able over the full redshift interval covered by the Lya for- 
est data (e.g., [Bouwens et al.[[2007l [Reddy fc Steidel[[2009 |. 
Where the data overlap, comparison of the Lya forest and 
the UV luminosity function allows us to constrain the es- 
cape fraction and limiting magnitude (minimum luminosity) 
down to which the luminosity function must be integrated 
in order to account for all the ionizing photons measured 
using the forest (Faucher-Giguere et al. 2008a). Since these 
parameters are constrained where the data overlap, we can 
test whether they must evolve with redshift in order to ac- 
commodate the reionization constraints from WMAP and 
galaxy surveys. Such evolution, in particular in the escape 
fraction, is sometimes invoked to support the hypothesis 
that galaxies can indeed reionize the Universe (e.g., Haardt[ 
fc Madau|2011 1, but there is little direct evidence for the re- 



quired change because direct measurements of escaping Ly- 
man continuum photons are prohobitive during the epoch 
of reionization. 

The plan of this paper is as follows. In ^ we review how 
UV luminosity function measurements can be converted into 
predictions for the reionization history. We show how uncer- 
tainties in the SED of galaxies, their escape fraction, and the 
limiting magnitude introduce large degeneracies and allow 
a wide range of scenarios to be consistent with the stan- 
dard WMAP constraint. In ^ we introduce the Lya forest 
constraints on the ionizing background at 2 z ^ 6 and 
explain how these constraints relate to the ionizing sources. 
In Q we compare with the galaxy UV luminosity function 
and Lya forest data to constrain the escape fraction and 
limiting magnitude at z — 4. We then combine these con- 
straints with the higher-redshift galaxy survey data and the 
measured WMAP optical depth to quantify the allowed sce- 
narios, parameterized by the required limiting magnitude 
and evolution of the escape fraction. We conclude with a 
discussion of the implications for galaxy formation and ex- 
periments aimed at probing the epoch of reionization in ^ 

Throughout, we assume cosmological parameters 
consistent with the WMAP 7-year data in combina- 
tion with supernovae and baryonic acoustic oscillations: 
(fim, rib, ^A, h) = (0.28, 0.04 6, 0.72, 0.7) (abbreviated 
WMAP-7; IKomatsu et al.[[2011[). We adopt hydrogen and 



helium mass fractions X = 0.75 and Y = 0.25, respectively. 
All magnitudes are in the AB system ( [Oke fc Gunn][1983[ ), 
Unless otherwise noted, all errors are la. 



^ Because reionization is predicted to be highly inhomogeneous, 
existing constraints have not ruled out that some regions of the 
Universe may have been reionized as late as 2 ~ 5 (e.g., [Mc-[ 
[Greer et al.[[2011^ , but as we show in this paper various data 
taken collectively suggest that the bulk of reionization occured 
significantly earlier. 
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Figure 1. Fits to Schechter galaxy UV luminosity function pa- 
rameters versus redshift. Data points with error bars are from 
[Bouwens et al.| | |2011c^ . The solid line is our best fit linear model 
(FIT). The dashed and dotted lines show the MIN and MAX 
models, in which the parameters were adjusted within the linear 
fit formal Icr errors to minimize (MIN) or maximize (MAX) the 
contribution from faint galaxies (see text for details and Table [l] 
for numerical values). 



2 GALAXY SURVEY AND WMAP 
REIONIZATION CONSTRAINTS 

The two most basic observational constraints on hydro- 
gen reionization are the high-redshift galaxy UV luminosity 
function (LF) and the Thomson (electron) scattering optical 
depth to the microwave background measured by WMAP-7, 
Te = 0.088 ± 0.015 (|Komatsu et al.||2011|). In the follow- 
ing, 



we review how these measurements can be combined to 
constrain parameters of the ionizing source population, in 
particular the limiting UV magnitude Afum and the escape 
fraction of ionizing photons /esc from star-forming galaxies. 
The procedure is based on calculating, for a given set of 
assumptions on the galaxy population, the predicted evolu- 
tion of the IGM ionized fraction versus redshift and evalu- 
ating the corresponding Tc. As we show in |2.3[ the escape 
fraction is degenerate with the ratio of 1, 500 A UV contin- 
uum to ionizing flux, a quantity sensitive to the SED of the 
galaxies and whose effect we encapsulate in a dimensionless 
parameter ^ion defined below. Thus, our analysis formally 
constrains the combination (^ion/csc. For simplicity, though, 
we will occasionally summarize our results in terms of /esc 



(for values of ("ion motivated by stellar population synthesis 
models), since it is the most uncertain of the two factors. 

In this paper, we assume that the majority of the ion- 
izing photons are produced by star-forming galaxies domi- 
nated by ordinary Pop II stars. In principle, other sources 
such as massive Pop III stars (e.g., Bromm et al. 1999 



Yoshida et al. 20041, accreting black holes (e.g., Haiman 



& Loeb 1998 Madau et al. 2004 



Kuhlen & Madau 



20051, 



Belikov & Hooper 



20091 



or annihilating dark matter (e.g. 
could also contribute ionizing photons. However, there is 
essentially no observational support for these more exotic 
scenarios. In particular, the luminosity function of luminous 
quasars drops sharply beyond z ~ 2 (Hopkins et al. 20071 



and theoretical models suggest that only one or two super- 
novae from Pop III stars suffice to trigger the transition to 
Pop II in an early halo ( Wis e et al.|2010 1 . This is supported 
by IGM metallicity measurements at 2 — 5 — 6, which show 
that the relative abundances are consistent with measure- 
ments down to 2 ~ 2, and thus that there is no evidence for 
significant metal production from Pop HI stars in the first 
billion years ( Becker et al.|2011b l. In contrast, star-forming 
galaxies are now routinely observed at 2 > 7 and we show 
explicitly in this work that scenarios in which they are solely 
responsible for hydrogen reionization are consistent with the 



available data (for a recent review, see also Robertson et al, 
|2010|). 

If sources other than star-forming galaxies dominated 
hydrogen reionization, then the constraints on and /esc 
that follow would be arbitrarily weakened. However, follow- 
ing Occam's razor, we do not consider such scenarios further 
here. 

2.1 Calculation of the HIT volume filling fraction 
and of the Thomson optical depth 

The evolution of the volume filling fraction of ionized hy- 
drogen, Qmi{z), is given by the differential equation 

^QhII _ T^ior ^ 



(1) 



dt TlH tree 

consisting of a source term proportional to the ionizing emis- 



sivity and a sink term due to recombinations ( Madau et al 



19981 



Under the assumption that galaxies provide the bulk of 
the ionizing photons, the comoving ionizing emissivity (in 
units of photons per unit time, per unit volume) can be 
expressed as an integral over the galaxy UV LF, (/)(Muv): 



dMuV '^(A/uv)7ion(Muv) /es 



(2) 



We denote by 7ion(Muv) the ionizing luminosity (in units of 
photons per unit time) of a galaxy with absolute rest-frame 
UV (1500 A) magnitude Afuv. /esc denotes the effective es- 
cape fraction, which by definition we treat as a function of 
z only (see |2.4[ ). The volume averaged recombination time 
is given by 

1 



CraiaB(To) nH(l + Y/4X) (1 + z) 



0.93 Gyr 



To 



2 X 104 K 



1 + 2 

7 



(3) 

-3 



where as is the case B hydrogen recombination coefficient. 
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To is the IGM temperature at mean density, Chii is the 
efTective clumping factor in ionized gas, and nu is the mean 
comoving hydrogen number density. We assume that hehum 
is singly ionized at the same time as hydrogen, but only fully 
ionized later through the action of quasars (e.g., Faucher- 
Giguere et al. 2008a). We use the effective clumping factor 
to account for both the actual dumpiness of the gas and for 
the fact that the IGM temperature (and hence the proper 
recombination coefflcient) in general depends on density, so 
that formally an average over the temperature distribution 
should be performed. 

The clumping factor must be selected with care, since 
formal averages Chii ~ {'T'Hii)/{"hii)^ over simulation vol- 
umes yield large values ~ 30 (e.g., Gnedin &: Ostriker|1997 
Springel fc Hernquist|[2003 1 that imply very demanding re- 
quirements on the ionizing sources. These large clumping 
factors arise because the average includes very dense galaxy 
halo gas. However, absorption of ionizing photons by gas in- 
side (or in the immediate vicinity of) galaxies is already ac- 
counted for by the escape fraction. Thus, the correct clump- 
ing factor to use is one that accounts only for recombinations 
occurring in the more diffuse IGM. Although some ambigu- 
ity is inherent in this definition, recent studies suggest that 
values Chii = 1—3 are appropriate during the epoch of reion- 
ization (e.g ., [Pawlik et al.|2009| |Shull et al.||2MT| |McQuinn| 
et al.|20TT l. The IGM temperature To is also uncertain, but 
the fiducial value To = 2 x 10* K is reasonable for freshly 
reionized gas ( Hui fc Haiman|2003 1. 

The Thomson optical depth to microwave background 
is then obtained by integrating Qhii, 



r°° cd + z)'^ 

/ dz Qnu{z)aTnH{l + r,Y/4X), 

Jo -"v^J 



(4) 



where H{z) is the Hubble parameter, ctt is the Thomson 
cross section, and we consider helium to be only singly ion- 
ized {rj = 1) at z > 4 and doubly ionized {rj = 2) at lower 
redshift. 

The main uncertainties in these calculations, which we 
discuss next, are 

(i) the extrapolation of the LF to magnitudes and red- 
shifts for which no direct measurement exists (§ 



2.21 



(ii) the conversion from Muv to ionizing photon luminos- 
ity (7ion; §[23|, 

(iii) and the escape fraction of ionizing photons (/esc; 



2.4 1. 



2.2 High redshift galaxy luminosity functions 

We base our analysis on recent observational determinations 
of the rest-fram e UV LF at .z ^ 4 in the HUDF09 (|Beck- 
with ct al. 20061 lOesch et al.||2007|), ERS (IWindhorst et al. 



2011D and CA NDELS fields ([Grogin et al.|2011||Koekemoer 
et al.|2011| by |Bouwens et al.| ( |2007[|2011b|c| ) The best-fit 
Schechter function parameters {(j)* , M* , and a) are summa- 
rized in Bouwens et al. (2011c hereafter Bll). M* quanti- 
fies the characteristic magnitude, 0* measures the comoving 
number density, and a is the faint-end slope. 

To interpolate between redshift bins and extrapolate to 
redshifts not directly probed by the data, we fit the redshift 
evolution of the three Schechter parameters to a simple lin- 
ear model of the form {M* , logjo <^*, a} = A + B{z — 6). 
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Figure 2. Comparison of our FIT, MIN, and MAX luminosity 
function models to the data of [Bouwens et ah] ||2011a[ l at z = 4 
and 7, and to the updated limits from [Oesch et aT] | |2011| l at 
z = 10. The 2 = 10 data points, obtained from a single galaxy 
candidate, were not included in the fits. 



The best-fit parameters, denoted FIT, are given in Table [T] 
In order to explore the uncertainties in the extrapolation 
to very faint galaxies, we adopt two additional models, in 
which we vary the redshift evolution within the formal la 
errors of our linear fits to either maximize (MAX) or mini- 
mize (MIN) the contribution from faint galaxies. Compared 
to our FIT model, the MAX model has a slightly brighter 
and less rapidly dimming M* , a slightly larger and more 
slowly decreasing (j)* , and a steeper and more quickly steep- 
ening faint end slope a; and vice-versa for the MIN model. 
Figure [T] shows our three fits for the redshift evolution of the 
Schechter parameters. 

Figure [2] shows how these fits compare to the actual 
LF data from [Bouwens et al. (201 la I at z = 4 and z — 7, 
and to the la upper limits at z ~ 10 obtained from the 



detection of a single galaxy candidate by Oesch et al. ( 2011 1. 
While the FIT and MIN models are in good agreement with 
the 2; ~ 10 limits, the MAX model predicts more galaxies 
than observed at Muv ~ —19.6 by a ~ 2.5a. Given the 
substantial uncertainties in estimating limits from a single 
candidate in a relatively small field, we however consider the 
MAX model to represent a valid limiting case. 

For comparison with theoretical predictions, it is useful 
to relate the UV magnitudes to the total mass of the haloes 
likely to host these galaxies. Since direct mass determina- 
tions from gravitational lensing or clustering are not avail- 
able at very high redshift, we attempt to establish such a 
relation via an abundance matching technique, by equating 
the cumulative dark matter halo mass function to the cumu- 
lative UV luminosity function over the redshifts of interest. 
At lower redshift the validity of the abundance matching 
technique has been demonstrated by its ability to reproduce 
the spatial clustering of galaxies in the SDSS/LRG catalog 
( jConroy et ar][2006l [Moster et aT][20Tol |Guo et al.|[20Tol ). 
Here, we use the UV luminosity, which traces star forma- 
tion rather than stellar mass, and we should expect a larger 
scatter in its relation to total halo mass. Nevertheless, the 
relation is likely to still be monotonic on average and we 
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Table 1. UV luminosity function evolution models 



Model 






logio 


0* 


a 






A 


B 


A 


B 


A 


B 


FIT 


-20.42 ±0.05 


0.27 ± 0.03 


-3.01 ± 0.04 


-0.07 ± 0.02 


-1.84 ±0.04 


-0.06 ±0.02 


MIN 


-20.37 


0.30 


-3.05 


-0.09 


-1.80 


-0.04 


MAX 


-20.47 


0.24 


-2.97 


-0.05 


-1.88 


-0.08 



FIT denotes the best linear fit of the form {M* , logj^Q <f>* , a} = A + B(z — 6) to the Schechter parameters reported in |Bouwens et aT] 
( |2011c^ at 2 = 4, 5, 6, 7 and 8. The parameters of the MAX and MIN models are adjusted within la of the best fit (independently) so 

as to maximize and minimize the contribution of faint galaxies. 




-16 -15 -14 -13 
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Figure 3. Abundance matching between the dark matter halo 
mass function and the UV luminosity functions from Bouwens et 
al. (2010a) at 2 = 4, 7 and 10. 



therefore expect the abundance matching results to be valid 
at the order-of-magnitude level. The results are shown in 
Figure. [3] which reveals that the faint values of Afiim advo- 
cated by Bll correspond to total halo masses below 10^ Mq. 
These results are in good agreement with a similar determi- 



nation by Trenti et al. ( 2010 1 



2.3 Conversion from UV magnitude to ionizing 
luminosity 

To evaluate equation ([2|, it is necessary to convert from the 
measured UV magnitudes to ionizing luminosity (the 7ion 
term). To do so, we adopt a simple double power-law model 
for the galaxy SED in the relevant range|^ 



hv <lKy 

1 Ry s; ft// < 4 Ry 
hv > A Ry. 

(5) 



^ Some authors first convert the UV magnitude to a star forma- 
tion rate, then convert the star formation rate to a rate of pro- 
duction of ionizing photons. This however introduces extraneous 
steps. 



This simple form is adequate to capture the main features of 
more detailed stellar population synthesis models over the 
limited energy range of interest (cf. [Leitherer et al.||1999| 
Schaerer 20031. We use the notation u\ to denote the fre- 



quency corresponding to wavelength A/A, e.g. ^'912 is the fre- 
quency at the 912 A Lyman edge. Note that the "/5" slopes 
often discussed in the literature (e.g., Bouwens et al.|2010b \ 
are usually defined in terms of wavelength, L\ oc v~^'' , so 
that we have the relation = + 2). 

The hydrogen ionizing photon luminosity (7ion) is then 
given by 



7ion = 



dv 
hv 



2xl0^^s" 



ergs~^ Hz" 



Cion. (6) 



To express this as a function of UV magnitude, we use 
the standard AB relation logj^Q(L^^5oo/(ergs~^ Hz~^)) — 
0.4(51.63 — Muv). Using equation |5] we can solve for the 
dimensionless parameter C,ioii'- 



Cion — 1.5 



■fLyC 
0.2 



(1.65) 



- 4" 
7 



(7) 



This parameter, a function of the stellar spectrum character- 
istics, encapsulates all the information necessary to convert 
from UV magnitude to ionizing photon luminosity. 

In order to bracket the uncertainties in the spectral pa- 
rameters (/LyC, /Sf, 7), we consider three different models: 
a fiducial model (FID) with Cion = 1, a harder spectrum 
model (HARD) with Cion = 2, and a softer spectrum model 
(SOFT) with Cion ~ 0.5. This range is representative of Pop 
II star-forming galaxies with continuous star formation his- 
tories and age ~ 10 - 100 Myr ( [Leitherer et al.|1999| ). Note 
that converting from Muv to ionizing luminosity via the 
star formation rate as done in Bll corresponds to Cion = 1 
(our FID model) . Our three Lv models thus span a factor of 
two variation (up and down) around the hydrogen-ionizing 
luminosity used by Bll. 



2.4 The escape fraction of ionizing photons 

Some fraction of the ionizing radiation produced by stel- 
lar populations is absorbed by dust and neutral hydrogen 
within their host galaxies, and thus does not contribute to 
ionizing the IGM. We capture this suppression by a simple 
multiplicative prefactor, /osc, applied in equation ([2|. Since 
our calculations are tied to the observed rest-frame UV LF, 
our /esc is strictly speaking a relative escape fraction, cap- 
turing the additional suppression of photons blueward of the 
Lyman edge compared to 1500 A photons. While neutral hy- 
drogen only absorbs the ionizing photons, dust extinguishes 



© 0000 RAS, MNRAS 000, 000-000 



6 Kuhlen & Faucher-Giguere 



LF 


L„ 


fee 


M,„ 




FIT 


FID 


0.20 


-10 


0.079 


FIT 


HARD 


0.20 


-14 


0.073 


MAX 


HARD 


0.20 


-16 


0.074 


MAX 


HARD 


0.05 


-12 


0.080 


... MIN 


HARD 


O..^ 


-14 


0.074 


MAX 


SOFT 


0.50 


-14 


0.084 


MAX 


FID 


0.50 


-16 


0.079 




14 16 18 



Figure 4. Volume fraction filling of HII regions as a function of 
redshift for a set of representative models that satisfy the mea- 
sured galaxy UV LF and the WMAP-7 Thomson scattering op- 
tical depth. LF evolution fits FIT, MIN, and MAX are shown 
in blue, cyan, and magenta, respectively. The FID, SOFT, and 
HARD spectral hardness models are indicated with solid, dotted, 
and dashed lines. The line thickness corresponds to /esc = 5%, 
20%, and 50% (from thin to thick). 



1500 A and ionizing photons similarly. Because of this broad 
band extinction by dust, /esc is not equal to the fraction of 
all ionizing photons produced by stars which are absorbed 
in the galaxy. Evaluating the latter would require knowl- 
edge of dust extinction, but is not actually required for our 
purposes. Similar relative definitions of the escape fraction 



are often adopted observationally as well (e.g., |Steidel et al 
2001[|Shapley et al.||2006| |Inoue et al]|2006[ ) 



The true escape fraction may well vary with galaxy 
mass, age, star formation history, or other properties. Such 
dependences are however essentially unknown at this time. 
We therefore assume in this work that fasc is a function of 
z only, i.e. we use /esc (2) to represent an effective escape 
fraction averaged over the galaxy population at redshift z, 
suitably weighted by the (unabsorbed) ionizing luminosity. 
A time dependence of /osc could thus arise from either a gen- 
uine time evolution in the escape fraction of galaxies (e.g., 
owing to an evolution in the star formation rate and its asso- 
ciated feedback), or from a redshift evolution in the make up 
of the galaxy population, with the escape fraction of galaxies 
with certain properties remaining constant. In !|4] we quan- 
tify the redshift evolution required of /esc required by the 
data for different scenarios. 



2.5 Range of models allowed by the UV LF and 
WMAP-7 constraints alone 

In Figure [5] we show the Thomson optical depth, Tc, for 
different reionization scenarios consistent with the measured 
UV luminosity function. The models explored correspond 
to varying assumptions for Mum, Cion, and /esc, which are 
further assumed here to be constant with redshift. 
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Figure 5. Thomson scattering optical depth to the microwave 
background versus limiting UV magnitude. The colors represent 
our three different galaxy UV LF parameterizations: FIT (blue), 
MIN (cyan), and MAX (magenta). The solid line corresponds to 
the FID (Cion = 1) Li/-model, and the shaded regions are bounded 
by the SOFT (Cion = 0.5) an d HARD (Cion = 2) models. The 
WMAP-7 Te = 0.088 ± 0.015 jKomatsu et al.|2011| is indicated 
with a gray band. The top panel is for /esc = 20%, the bottom 
left for /esc = 5% and the bottom right for /esc = 50%. M^^, 
Cion I a-iid /esc are assumed constant in these calculations, and we 
used a clumping factor of Chii = 3. 



For the best-fit UV LF evolution parameterization 
(FIT), fiducial L^-model (Cion = 1, FID), and /osc = 0.2 
(solid blue line in the top panel), we recover the result of 
Bll that a very faint limiting magnitude, Miim ^ — H, is 
required in order to produce an optical depth in agreement 
with WMAP-7. However, many other solutions are possible. 
For example, the same LF model with a harder spectrum 
(upper edge of blue shaded region) is consistent with the 
WMAP-7 data for Mum = -14, and with the MAX LF 
model (magenta band) the WMAP-7 Tc constraint can ac- 
commodate values of A/um ranging from —11 to as bright as 
— 16, depending on the spectral hardness. The escape frac- 
tion provides yet another degree of freedom. With a constant 
/esc of 5% (bottom left panel), most models cannot satisfy 
the WMAP-7 constraint. On the other hand, if /esc = 50% 



© 0000 RAS, MNRAS 000, 000-000 



Concordance models of reionization 7 



even the MIN LF fit or models with very soft spectra can 
result in a sufficiently high Tc. 

In Figure |4] we show the volume filling factor of HIT 
regions, Qmi{z), for a few representative models that all 
satisfy the WMAP-7 Tc constraint. Interestingly, these sce- 
narios have limiting magnitudes ranging from —10 to —16 
and include models that extrapolate the contribution of 
faint galaxies quite differently. We conclude that the existing 
measurements of the high-redshift galaxy luminosity func- 
tion (still limited to relatively luminous sources) and of the 
Thomson optical depth to the microwave background do not 
uniquely determine how reionization proceeded. In particu- 
lar, these constraints do not suffice to determine the role 
played by low-luminosity galaxies. 



This mean transmission measurement, based on 86 high- 
resolution and high-signal-to-noise quasar spectra covering 
LyQ redshifts 2 ^ z ^ 4.2, was corrected for absorption by 
metal ions and for biases in the continuum fits, an important 
effect at 2 > 4. At z = 5 and z = 6, we use the constraints 



on Phi from Bolton & Haehnelt (20071, also from mean 



transmission data. We do not use proximity effect measure- 
ments, as they are typically of lower statistical precision and 
affected by more severe systematics (e.g., Faucher-Giguere 
et al. 2008c). Nevertheless, at z = 5 — 6, where some of 
these effects are mitigated, the proximity effect measure- 



ments of Calvcrlcy et al. (20111 are consistent with Bolton 



fc Haehnelt, (2007) . 

In Table |2| we summarize the Phi measurements and 
other inputs used in our Lya forest analysis. 



3 Lya FOREST CONSTRAINTS ON THE 
IONIZING SOURCES 

The Lya forest provides complementary constraints on 
the cosmological emissivity of ionizing p hotons ([Miralda- 
Escude"2003] [Bolton fc Haehnelt] |2007| |Faucher-Giguere 
et al.||2008ai ). Although saturation prevents accurate mea- 
surements of the Lya forest at z > 6 (e.g.. Fan et al.|2002 1, 
it has the advantage of being a complete probe, in the sense 
that it includes the contribution of all ionizing sources, even 
if they are individually too faint to be detected in galaxy 
surveys. Furthermore, the ionizing emissivity implied by the 
Lya forest does not depend on an assumed escape fraction. 
Thus, the Lya forest is not subject to the two main un- 
certainties affecting the inference of the ionizing emissivity 
from the galaxy UV luminosity function, namely Mum and 
/esc. By assuming continuity between the post-reionization 
epochs probed by the Lya forest and the reionization epoch 
probed by high-redshift galaxy surveys, it is therefore pos- 
sible to significantly reduce the permitted parameter space. 
In particular, comparison of the UV luminosity function and 
Lya data where they overlap allows us to constrain a com- 
bination of escape fraction, the limiting magnitude, and the 
conversion factor from 1,500 A UV to ionizing luminosity 



3.1 Total ionization rate from the Lya forest 

The basic quantity constrained by the Lya forest is the hy- 
drogen photoionization rate Phi, 



Phi(z) = 47r 



du 
hv 



(8) 



where is the average specific intensity of the ultra-violet 
background, omiy) is the photoionization cross section of 
hydrogen, and the integral is from the Lyman limit to infin- 
ity. Indeed, the mean level of transmission of the Lya forest 
is set by the equilibrium between the ionizing background 
and recombinations in the IGM. Thus, given a model of the 
density fluctuations in the IGM and knowledge of the inter- 
galactic gas temperature-density relation, the mean trans- 
mission of the Lya forest can be inverted to give Phi (e.g., 
Ranch et al.|1997 l. 

In this work, we use principally the Phi data points from 
Faucher-Giguere et al. (2008a,b) based on the mean trans- 
mission measurement of Faucher-Giguere et al. (2008d). 



3.2 Prom ionization rate to ionizing emissivity 

The quantity most directly related to the sources of ioniz- 
ing photons is their spatially-averaged emissivity, (here 
in proper, specific units). Assuming that the ionizing back- 
ground has a power-law spectrum Ji, = J^g^a (^^71^912) "''^ 
between the HI and Hell ionizing edges (and zero beyond), 

1 

1 - . (9) 



PHi/t(7bg + 3) 
47rcrHi (1^912) 



1 

4?n^^ 



Since 



47r 



JAz) 



Amfp(«^, 2)' 



(10) 



where A^fp is the mean free path of ionizing photons in 
proper units (denoted A^fp at the Lyman limit) Q 

PHi(z)ft(7bg 



-3) 



'-HiAii^( 



4'>'bg 



+3 



(11) 



Assuming similarly that e„ = ei/912 ('^/t'9i2) ~' between the 
HI and Hell ionizing edges (and zero beyond)]^ 

du 



\z)^ 



1 



{l + zf 

1 



(12) 



'^912 

Vm{z)h 



(7bg + 3) 



{1 + zf ani>?^l{z) 



1 - 



1 

47 



47bg+3 

(compare with eq. ([2|). Equation (12 1 shows how the total 
comoving emissivity of ionizing photons can be derived from 
the photoionization rate measured from the Lya forest and 
knowledge of the mean free path of the ionizing photons, 
given a model for the spectral shape of the ionizing sources 
and their integrated background. 

At 3.6 ^ z ^ 4.2, we use the mean free path measured 



by Prochaska et al. (20091 using a stacking analysis. This 



approach avoids the usual uncertainties in calculating the 
mean free path from the column density distribution stem- 
ming from the difficulty of measuring the column density of 



This approximation is valid at 2 ^ 2, where the mean free path 
is much smaller than the Hubble scale. 

^ Since the ionizing background spectrum is affected by filtering 
by the IGM, in general 7 ^ 7bg (e.g., [Haardt fc Madau|[l996[ 
Faucher-Giguere et al.|2009 i. 
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Table 2. Lya forest constraints on the ionizing emissivity 



z 


Thi 


A912 

mrp 


^com 
"ion 


References 




10-" s-i 


pMpc 


10^0 s-1 cMpc-3 




9 fl 






Z.UltU.O l,_i,4J 




9 9 


n f^i -1-0 1 n 


997-1- (^1 






9 4 


u.oumu.uo 


1 ( '4:111^0 






9 

Z.D 


n f^i -1-0 (17 




9 0-1-0 f'^-'--^^ 
Z.UmU.U 




9 Q 


n cri _|_n fif:; 


iUDitzu 






3.0 


n 5Q-|-n 07 


84.4di21 






3.2 


0.66±0.08 


67.9±16 


3.3±0.9 (tl'^) 


FG08, scio 


3.4 


0.53±0.05 


55.2±13 


2.8±0.7 


FGO8, scio 


3.6 


0.49±0.05 


49.5±2.1 


2.6±0.3 


FGO8, P09 


3.8 


0.51±0.04 


41.7±2.4 


2.8±0.3 (t^;^) 


FG08, P09 


4.0 


0.55±0.05 


34.0±3.1 


3.2±0.4 itll) 


FG08, P09 


4.2 


0.52±0.08 


26.2±3.9 


3.5±0.8 (t^;^) 


FG08, P09 


5.0 


n 59+0 35 


13.9±3.6 


4.3±2.6 (±2.6) 


B07, SCIO 


6.0 


<0.19 


7.0±2.0 


<2.6 (< 2.6) 


B07, SCIO 



The HI photoionization rates measurements are taken from Faucher-Gigere et al. (2008a) (FG08) an d [Bolton &: Haehnelt] | |2007| (B07); 

the mean free paths are taken from the fits of |Prochaska et al.| l |200^ (P09) and |SongaiIa fc Cowie| l |2010[ l (SCIO). Errors on Fhi and 
"^mfp ^'^^ predominantly statistical (except for the B07 Fhi points, which include a systematic error budget). Total uncertainties 

on 5^^°™, including systematic effects arising from the spectral shape of the UV background and the thermal history of the IGM, are 
given in parentheses and shown by the light gray band in Figure |6] (see the text). The prefixes 'p' and 'c' indicate proper and comoving 

units, respectively. 



systems near the Lyman limit (on the flat part of the curve 
of growth). At the other redshifts 2 ^ 2 ^ 6, we use the 



mean free path derived by Songaila & Cowie (20101 based 



on a new analysis of the column density distribution. These 
expressions agree well with the mean free path inferred pre- 
viously by Faucher-Giguere et al. (2008a), but have signifi- 
cantly reduced uncertainties. On the other hand, this mean 



free path is larger than that assumed by Madau et al. ( 1999 1 
by a factor ~ 2.5. Furthermore, these mean free path mea- 
surements are significantly more accurate than the simple 
model based on the mean spacing between Lyman limit sys- 



tems assumed by Bolton & Haehnelt ( 2007 \ 



Figure [6] summarizes the IGM observational constraints 
on nfo™. The error bars on the data points account for the 
statistical uncertainty on the photonization rate and on the 
mean free path. The z — 5 and z = 6 error bars also include 
a systematic error budget on Phi, as quantified by |Bolton| 
& Haehnelt (20071. Total uncertainties, including systemat- 



ics, are indicated by the Ught gray band and estimated as 
follows. First, we allow for a 50% systematic error on the 
Fhi data points from Faucher-Giguere et al. (2008a) to ac- 
count for uncertainties in the thermal state of the IGM and 
the probability distribution function of density fluctuations, 
which enter in the mean transmission method (Bolton et 
al. 2005; Faucher-Giguere et al. 2008a). Second, we vary the 
source spectral index 7 from 1 to 3. The harder value 7 = 1 is 
preferred by optical line ratio diagnostics in local starbursts 
| |Kewley et al.|2001| ), while many stellar population synthe- 



sis models predict 7 « 3 (e.g., Leitherer et al. 19991. This 
range of slopes is also consistent with the possibility that 
quasars, with mean spectral index ~ 1.6 ( [Teller et al.|2002[ ), 
contribute significantly at the lower redshift end. The fidu- 
cial value assumed in our calculations is 7 = 1; since stellar 
population synthesis models generally predict softer spectra, 
we do not explore harder values. The spectral index of the 



background in the ionizing regime is not independent but 
instead satisfies 7bg = 7 — 3(/3 — l)[^where /3 is the slope of 
the HI column density distribution (Faucher-Giguere et al. 
2008a). We adopt /3 = 1.3 ( [Songaila fc Cowie|2010l ). 

To indicate the total uncertainty, we first calculate the 
range of nfo™ values allowed by simultaneously varying the 
systematically uncertain parameters to their extremes. We 
then add the statistical uncertainty to the minimum and 
maximum values in each redshift bin. We believe that this 
procedure conservatively captures the constraints on nf°™. 

The mean transmission data points are instantaneous 
constraints that must be satisfied by the galaxy popula- 
tion. Models of the ionizing background at intermediate red- 
shifts indicate that star-forming galaxies dominate at z > 3 
(Haehnelt et al. 2001; Bolton et al. 2005; Faucher-Giguere 
et al. 2008a, 2009). However, quasars may dominate the hy- 
drogen photoionization rate at later times, so that the total 
background should be regarded as an upper limit to the 
contribution of star-forming galaxies alone. It should also 
be noted that any redshift evolution in nf°™ contained in 
the light gray band in Figure [6] is allowed. In particular, it 
is possible that the true redshift evolution is titled in slope 
relative to that suggested by the fiducial data points. This is 
because the uncertain parameters could evolve significantly 
with redshift. For instance, measurements indicate that the 
IGM temperature peaks at To > 2 x 10'' K around z ~ 3.4 



^ This approximation (arising from the frequency dependence of 
the mean free path; eq. \W\ ) is valid at least up to 2 ~ 4, where 
the column density distribution has been measured to be well 
approximated by a series of power laws jProchaska et al.|2010| . 
However, it may break down at earlier times, where optically thick 
absorbers could be relatively more numerous and dominate the 
opacity. This provides additional motivation for focusing our post- 
reionization analysis at 2 = 4 (iQ. 
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Figure 6. Lya constraints on the rate at which ionizing photons are injected into the IGM (see jjsjfor details and Table[2]for numerical 
values). The light gray band indicates instantaneous constraints from the measured mean transmission of the Lya forest, including 
systematic effects. The dark gray band indicates the minimum value necessary to keep the Universe ionized, assuming that reionization 
is complete, for a fiducial IGM temperature Tq = 2 X 10* K and effective clumping factor Chii = 1 — 3. Models of the ionizing background 
indicate that the ionizing emissivity is dominated by star-forming galaxies at z > 3, but that quasars may dominate at lower redshifts 
(Faucher-Giguere et al. 2008a, 2009). 



(Lidz et al. 20101, possibly owing to re-heating from Hell 
reionization, but could be less than ~ 10'* K between HI 
and Hell reionizatio n (|Hui fc Haiman|[2003j [Becker et al 



2011a Bolton et al. 2012 1. Recent observations also indi- 
cate that the UV slopes of 2 ~ 7 galaxies are significantly 
bluer than their 2 ~ 3 counterparts | |Bouwens et al.|2010b[ ), 
so that the relevant spectral indexes could also evolve. 

It is apparent from Figure [6] that the best-fit comov- 
ing ionizing photon emissivity, nf"^, increases from 2 = 2 
to 2 = 4.2, and perhaps even to higher redshift. This basic 
behavior stems from the fact that while the photoioniza- 
tion rate, Fhi, is approximately constant over this redshift 
interval, the mean free path decreases rapidly with increas- 



ing redshift (Faucher-Giguere et al. 2008a McQuinn et al 



2011 1. Thus, an increasing ionizing emissivity is required to 
maintain the observed photoionization rate (eq. |2|. 



3.3 Keeping the Universe ionized 

Figure [g] also shows (as the dark gray band) n^°^''^^'^, the 
minimum nfo ™ required to keep the IGM ionized once it has 
already been reionized. This number is obtained by balanc- 
ing the global recombination rate in the fully ionized IGM 
with the rate at which ionizing photons escape galaxies: 



CHiiaA(ro)nH(l + y/4X)(l + 2) 
cMpc- 



(13) 



To 



2 X 10* K 



-0.7 / , , X 3 

1 + 2 



Here, oa is the case A recombination coefficient of hydrogen. 
Although we used the case B coefficient for the calculation 
of the HII volume filling factor during reionization (eq. (|4|), 
a large fraction of the recombinations directly to the ground 
state at later times do not actually contribute to the ionizing 
background ( [Faucher-Giguere et al.|2009| ). 



The nl°^''^^^^ band in Figure [H] covers the range Chi = 
1 — 3 and assumes a fiducial IGM temperature To = 2 x 10* 
K. This temperature is consistent with measurements at 
(e.i 



Lidz et 



al. 2010 1, but could be higher 



2 ~ 2 - 4 
or lower by 

re-ionizing sources and the time elapsed since reionization 



10 K depending on the spectrum of the 



(Hui & Haiman 20031. Estimates of the IGM temperature 
at 2 > 4 in fact suggest To < 10, 000 K at 2 ~ 5 - 6 ( [Becker 
[et al.[[201l"a| [Bolton et aL|[2012[ ). Fortunately, this uncer- 
tainty does not substantially affect our arguments, since 
the recombination coefficient is a relatively weak function 
of temperature. 



Note that the n^°^''^" lower limit only applies after 
reionization is complete. Since we do not a priori know the 
redshift of reionization, the plotted constraint (extending 
to 2 = 7.5) need not necessarily be satisifed everywhere. 
However, any viable reionization scenario must satisfy this 
constraint at all redshifts following the time when an ionized 
fraction ~ 1 is reached. 
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Figure 7. Top; Value of Cion/csc at z = 4 needed to simulta- 
neously match the total comoving emissivity of ionizing photons 
measured from the Lyo forest, and the observed UV lumi- 

nosity function at the same redshift, as a function of the limiting 
UV magnitude. Since our LF fits are almost identical at z = 4, we 
only show the FIT case. Because ^jon = 1 for the fiducial spectral 
model, the values directly quantify the implied escape fraction. 
Bottom: Power-law index k of the redshift evolution of Cion/csc 
(see eq. |15[ l needed to simultaneously match the z = 4 Lya forest 
and WMAP-7 Thomson optical depth constraints, as a function 
of the limiting UV magnitude (assumed constant here), for our 
three LF evolution fits. The solid (dotted) lines correspond to 
the median value (itlcr) of the WMAP-7 Thomson optical depth. 
The shaded regions encompass the total (including systematic) 
uncertainty in »^f°™(2 = 4). Note that some models with bright 
Miinj do not admit solutions for the entire range. 



4 COMBINING LYa FOREST, GALAXY 
SURVEY, AND WMAP CONSTRAINTS 

We showed in i|2]that many different reionization scenarios 
are consistent with the existing galaxy survey and WMAP 
constraints, principally due to uncertainties in Mum and 
/esc, and also ^ion- We now combine these constraints with 
the Lya forest data at lower redshifts (*j3|, which allow us 
to break certain degeneracies and quantify possible redshift 
evolution in the relevant parameters. The key idea is that 
for any Afum, comparison of the Lya forest and galaxy UV 
luminosity function data where they overlap imply a unique 




0.03 



Figure 8. Evolution of Cion/csc versus redshift required to simul- 
taneously satisfy the z = 4 Lya forest and WMAP-7 Thomson 
optical depth constraints, for M\\-^ = —10, —13, and —16 (models 
corresponding to the solid lines in Fig.[7]|. A ceiling of Cion/csc < 2 
is imposed in our calculations, corresponding to /esc = 1 for 
Cion = 2 (our HARD spectral model). 



Cion/osc value (assuming that galaxies dominate the ioniz- 
ing background). We choose to make this comparison at 
z = 4, because at this redshift we expect the hydrogen ion- 
izing background to in fact be dominated by star-forming 
galaxies (Faucher-Giguere et al. 2008a, 2009). Furthermore, 
at this redshift the observational constraints on the Lya 
forest transmission and the mean free path of ionizing pho- 
tons are quite good. Additionally, we do not expect this 
redshift to be strongly affected by large inhomogeneities 
in the ionizing background. At higher redshifts (especially 
at z ~ 6), interpretation of the Lya forest data becomes 
more uncertain because of the small number of sight lines 
available and because reionization may not be 100% com- 
plete ( McGreer et al.|[2011 l. Because only certain combina- 
tion (Miim, Cion/csc) are allowed at z = 4, we can quantify 
whether redshift evolution in these parameters is required 
in order to simultaneously satisfy the higher-redshift con- 
straints from galaxy surveys and WMAP. 

In principle, Afum, Cion, and /esc can all be arbitrary 
functions of redshift. However, the limited data do not al- 
low us to explore the entire range of possibilities. Instead, we 
adopt simple, single-parameter power-law models to quan- 
tify the required redshift evolution: 



(Cion /esc ) (2') (Cion/csc) (2- 4) 



1 + 2 

5 



10 



-0.4A/,i,„(z) 



= 10 



-0.4Af,i„{2=l) 



1 



(14) 



(15) 



The last expression encodes a power-law parameterization 
in the UV luminosity, which is logarithmically related to 
the UV magnitude. In order to prevent unphysical values of 
/esc, we impose Cion/csc 5; 2, Corresponding to a ceiling of 
/esc = 1 for the case Cion = 2 (our HARD spectral model; 
[2.3 1. There is also a physical lower limit on the luminosity 



of the faintest galaxies, corresponding to the minimum halo 
mass in which baryons can collapse and form stars. However, 
this limit is not as well understood and could depend signif- 
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Figure 9. Redshift at which reionization is 20% (top left), 50% (top right), 90% (bottom left), and 100% complete (bottom right), as 
a function of the limiting UV magnitude, for models in which Cion/csc(^) is tuned to reproduce both the WMAP-7 Thomson optical 
depth and the 2 = 4 Lya forest constraints. The colors represent our three LF evolution fits, and the shaded region encompasses both 
the WMAP-7 re Icr region and the total (including systematic) uncertainty in nf°^{z = 4). 



icantly on redshift. We impose an extremely faint ceiling of 
Miim = (~ 4 orders of magnitude below the UV suppres- 
sion sacle of M = —10), and discuss the viability of different 
scenarios later. Note that we have defined the signs of the k 
and A power-law indexes so that positive values correspond 
to increasing efficiency of ionizing photon production going 
to higher redshifts. 



In the following, we consider possible evolution first in 
Cion/esc ( ^4.1[ ), aud theu in Miim ( ^4.2[ ). In principle, there 
could be simultaneous evolution in both (^ion/csc and Afum, 
but the data do not allow us to discriminate between such 
mixed scenarios. Furthermore, we will show that the strong 
redshift evolution required in models relying only on rel- 
atively bright galaxies is most plausibly accounted for by 
evolution in the escape fraction. In §4.3| we consider the 
redshift and duration of reionization in different allowed sce- 
narios, showing that the redshift of 50% ionized fraction is 
consistently at 2:reion(50%) ~ 10 among the different models 
allowed by the data, but that the duration of reionization, 

A^rcion = Zrcion(100%) — ZrGion(20%) (where 2:roion(a;) is the 

redshift such that QHii(zrcion(a;)) — x), is on the other hand 
sensitive to the contribution of faint galaxies. 



4.1 Redshift evolution of C^ion/csc 

We first focus on the redshift evolution of Cion.fcsc from z — 4 
toward higher redshifts. 

Our goal is to determine as a function of Afum (here 
assumed to be independent of z) what values of k are con- 
sistent with both the z ^ 4 Lya forest and WMAP-7 Thom- 
son optical depth constraints, while also satisfying the mea- 
surements of the galaxy UV LF at the bright end. Oper- 
ationally, we first determine for a given Mum what range 
of Cion/esc (« = 4) is required to give nf°^{z = 4) = 



3.2 



+2.2 



X 10^ 



cMpc" 



The results of this calculation 



are shown in the top panel of Figure [7] For A/um — —10 
to —16, the Lya constraints require i^ionfcsc{z — 4) to lie 
between 2% and 8% (since the luminosity function has al- 
ready been measured down to Afuv = — 16 at z = 4, cases 
with brighter A^um at this redshift are not allowed; IBouwens] 
et al.|2007 l. As Cion ~ 1 for the fiducial spectral model, this 
directly quantifies the implied escape fraction. This result is 
independent of which LF fit we employ, since at z — 4 the 
parameters of our three fits are nearly identical. 

In the second step, we determine for each Mum and 
Cion/csc(2 = 4) what valucs of the power-law index k yield a 
Thomson optical depth in the range allowed by WMAP-7. 
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Figure 10. Same as Figure|6] but with the models from Figure|4] 
overplotted, for the constant Cion/csc case (top) and modified to 
allow for redshift evolution in Cion/osc (bottom). For the models 
with redshift evolution in Cion/esc! the luminosity function pa- 
rameterization and Miim are fixed to the values from Figure 4, 
Cion/csc(2; = 4) is set by the Lya forest data at z = 4, and we solve 
for K (eq. such that the Thomson optical depth matches the 

central WMAP-7 measurement. In the same order as the legend in 
Figure|4] the best-fit values are k = 2.2, 3.3, 3.0, 1.3, 6.7, 2.1, 3.0. 
The model lines are greyed-out at z < 4, since we do not utilize 
galaxy luminosity function data at those redshifts. 



The resulting range of allowed k values is shown in the bot- 
tom panel of Figure [7[ with the different color bands corre- 
sponding to our three LF fits. The width of the bands encom- 
passes both the 1(7 uncertainty of the WMAP-7 Thomson 
optical depth measurement and the total (including system- 
atic) uncertainty in the z — 4 Lya forest data. Figure [s] 
shows the curves of (^ion/csc versus z corresponding to al- 
lowed values of k, for representative choices of i\fiim. For 
each K solution, Figure [o] shows 2;rcion(20%), 2;roion(50%), 

2reion(90%), and «rcion ( 100%) . 

Lastly, Figure [lO] shows explicit examples of how red- 
shift evolution in ("ion/osc allows models to simultaneously 
satisfy reionization-epoch constraints and the z < 6 Lya for- 
est data. Specifically, we consider the same models as in Fig- 
ure|4](top) and modify them to allow for redshift evolution in 
Cion/esc (bottom). For the models with redshift evolution in 
Cion/esc, the lumiuosity function parameterization and Mm-n 
are fixed to the values from Figure 4, Cion/esc (^ = 4) is set by 
the Lya forest data at z — 4, and we solve for k such that 
the Thomson optical depth matches the central WMAP-7 
measurement. While the original models without evolution 
in Cion/esc overproducc the ionizing emissivity probed by the 
LyQ forest at « = 4, the modified models simultaneously sat- 



isfy all the constraints. Furthermore, models exist in which 
the extrapolation to z = 2 is also in good agreement with 
the lower-redshift Lya forest data. 

At the 1(7 level, only the MAX model with very faint 
Afiini > —11 can accommodate no redshift evolution in 
Cion/osc. However, these scenarios are disfavored by exter- 
nal constraints on the duration of reionization from the ki- 
netic Sunyaev-Zeldovich effect (Zahn et al. 2011), which in 
combination with the WMAP-7 optical depth constrains the 
timing of its beginning and end. For the most conservative 
case of arbitrary correlations between the thermal Sunyaev- 
Zeldovich (tSZ) effect and the cosmic infrared background 
(CIB),|ZaMetaI]([201l]) find that 2reion(20%) < 13.1 at 95% 
confidence level (CL), and 2reion(99%) > 5.8, also at 95% 
CL. In this work, we take these constraints at face value. It 
is important to bear in mind, however, that the templates 
on which they are based assume that reionization occurs 
primarily via star-forming galaxies. Furthermore, the limits 
on the kSZ signal rely critically on accurate subtraction of 
contaminating point sources. It will thus be important to 
confirm these findings with refined analyses. 

For the best-fit parameterization of the UV LF (the 
FIT model), models with no redshift evolution in Cion/esc are 
disfavored even for — —10. Models that rely only on 

brighter galaxies formally satisfy all the present constraints 
but only for strong redshift evolution in Cion/esc- For exam- 
ple, the case of Afum = — 16 for the FIT parameterization 
requires an evolution in Cion/esc by a factor ~ 20 from z — 4 
to z = 9. As we will discuss at greater length in [|5j mod- 
els that rely too heavily on fainter galaxies may be in ten- 
sion with theoretical models that suppress star formation 



in early, low-mass systems (e.g., Krumholz & Dekel 2011 



Kuhlen et al. 20121, which are helpful in explaining some 



properties of the cosmic star formation history. If star for- 
mation is indeed suppressed in those early dwarfs, then the 
existing data would imply strong evolution in the escape 
fraction. 

We also explored constraints on the redshift evolution 
Cion/esc from z = 2 to z = 4 by comparing the Lya forest 
data to the galaxy UV luminosity function from |Reddy Ei\ 
Steidel (20091 a.t z — 2. Over that redshift interval, a wide 
range k ~ — 4.5 is allowed, almost independent of the 
assumed Mum owing to the relatively shallow faint-end slope 
of the luminosity function. In particular, the combination of 
the UV luminosity function and Lya forest data alone do not 
require any significant evolution. Note that such evolution 
is nonetheless allowed by the data, and in fact suggested 
by direct Lyman continuum observations (e.g., [lnoue et alT] 
2006 Siana et al.||2010 l. Although this is not necessary on 
physical grounds, it is interesting that most of the k values 
implied from z — 4 and up (Fig. [7| are also allowed from 
z = 2 to 2 = 4. 



4.2 Redshift evolution of Mi^-n 

We now turn to the possibility of redshift evolution in Miim. 
Starting from the same values of Cion/esc (z = 4) for a given 
Muin{z ~ 4) as in the previous section, we determine what 
values of A can produce agreement between the Lya forest 
constraints a,t z — 4, and higher-redshift LF and WMAP-7 
constraints. The top panel of Figure [ll] demonstrates that 
substantial evolution is necessary to match the WMAP-7 
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Figure 11. Top: Power-law index A of the rcdshift evolution of 
the limiting UV magnitude M\\^ (see eq. 14) needed to simulta- 
neously match the z = 4 Lya forest and the WMAP-7 Thomson 
optical depth constraints, as a function of M\\-^{^z = 4) (Cion /esc 
is assumed constant). As before, the solid (dotted) lines corre- 
spond to the median value (ztlcr) of the WMAP-7 Thomson op- 
tical depth, and the shaded regions encompass the total (includ- 
ing systematic) uncertainty in '^f°™(z = 4). Models with the MIN 
LF fit do not admit any solutions with evolution only in Mum. 
Middle: Redshift at which reionization is completed (Qhii = !)■ 
Most of the Miini-only evolution models do not complete reion- 
ization by z = 5.8, as required by the combination of kSZ and 
WMAP-7 data l |Zahn et aL||2011| l. Bottom: Median M^^ red- 
shift evolution (corresponding to the solid lines in the top panel) , 
for A/iijn = —10, —13, and — 16. Very faint limiting magnitudes 
A^lim 2> — 10 are also likely excluded on physical grounds. 



Thomson optical depth measurement. Many models actually 
do not allow for a solution: for the MIN LF fit, the contribu- 
tion of dwarfs is suppressed to such a degree that no amount 
of Miim evolution is able to raise To into the range allowed 
by WMAP-7. The FIT model only has solutions with very 
steep Miim evolution, requiring Mum significantly below —10 
at high redshift. Only the MAX models are able to provide 
solutions with a moderate amount of Mum evolution. Most 
importantly, as the middle panel of Figure [TT] demonstrates, 
the majority of the models with evolving A/um (but constant 
Cion/csc) do not reach complete reionization by z = 5.8, as 
required by the combmation of the kSZ and WMAP con- 
straints ( Zahn et al.|2011 i. 

Except for small, extreme corners of parameter space, 
it is therefore not possible to simultaneously match z = 4 
LyQ forest and reionization constraints, and complete reion- 
ization in time, by allowing only evolution in M\\m- This 
provides further evidence for the need for a significant red- 
shift evolution in ^ion/csc from z = 4 toward higher redshifts, 
as discussed in the previous section. In contrast, the data do 
not provide conclusive evidence for a significant redshift evo- 
lution in Miim (although some evolution may be expected on 
physical grounds). We will therefore concentrate the follow- 
ing discussion on scenarios with constant but evolving 
Cion/osc. It is important to keep in mind, however, that evo- 
lution in Miim could reduce the required amount of evolution 



4.3 Redshift and duration of reionization 

Interestingly, Figure [9] shows that the redshift of 50% ion- 
ized fraction is consistently at 2 ~ 10 between the dif- 
ferent allowed scenarios. This is consistent with the red- 
shift of instantaneous reionization implied by WMAP-7, 
Zrcion = 10.4 ± 1.2. However, the predicted duration of the 
reionization process is more extended for scenarios that in- 
clude a larger contribution from faint galaxies. This is a 
consequence of the shape of the LF: the brighter and closer 
to M* (the knee of the LF) Miim is, the fewer galaxies con- 
tribute to reionization at high z. The reionization process 
then quickly sets in once the exponential cutoff of the LF 
has shifted to bright enough galaxies that Mum galaxies be- 
come common. In contrast, with a faint Miim abundant faint 
galaxies contribute to reionization even at very high red- 
shifts, and therefore the overall evolution of the process is 
slowed down. 

As shown in Figure|9j the kSZ data favor relatively short 
reionization histories and thus disfavor models that rely too 
heavily on high-redshift dwarfs, as in the MAX models. A 
range of MIN and FIT models are however allowed. In par- 
ticular, our conclusions are consistent with theoretical anal- 
yses anchored in the predicted dark matter halo mass func- 



tion, rather than to the observed luminosity function ( Trenti 
et aL|2010[|Mesinger et al.|2011||Ciardi et al.|2011| ), indicat- 
ing that as of yet undetected galaxies must contribute signif- 
icantly to reionization. However, our analysis suggest alter- 
native possibilities when strong redshift evolution in i^ion/csc 
is allowed. 

Recently, two significant observational advances in us- 
ing astrophysical sources to probe the epoch of reionization 
have been reported. First, [Mortlock et al.| ( [2011[ ) discovered 
a luminous quasar at z = 7.085 in which the Lya trans- 
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mission profile is consistent with an IGM neutral fraction 
~ 10% at that time. Although this is compatible with most 
of the scenarios allowed by our analysis (Fig. [9| , excluding 
those relying on a maximal contribution from dwarf galaxies, 
alternative interpretations exist in the context of inhomoge- 
neous reionization ( Bolton et al.|[20TT K 

Second, recent surveys for Lya emitting galaxies at 
z ^ 6 have found a decreasing fraction of Lyman break- 
selected galaxies with detected Lya emissio n of rest-frame 
equivalent width ^ 20 A from z ~ 6 to z ~ 7 ( Schenker et al. 
2011[|Pentericci et al.|2011| |Ono et al.|2012[ ). The decline is 



such that existing models of Lya propagation through galac- 
tic winds and the intervening IGM indicate Qhii 50% at 
z ~ 7 ( Dijkstra et al.|2011 1. Such a high neutral fraction at 
2 = 7 is in tension with the concordance scenarios summa- 
rized in Figure [9] Furthermore, if reionization is essentially 
complete by z ~ 6 as the combination of kSZ and WMAP-7 
measurements indicate, and as suggested by Gunn-Peterson 
troughs in the Lya forest (e.g., [Fan et al.|2002[ ), then a very 
rapid evolution in the neutral fraction would be implied. 
These findings are not easy to reconcile, and highlight the 
need for more detailed modeling of Lya radiative transfer 
in order to fully exploit the new and upcoming high-quality 
data. 



5 SUMMARY AND DISCUSSION 

Measurements of the galaxy UV (~ 1,500 A) luminosity func- 
tion at redshifts z > 6, recently improved by more than an 
order of magnitude thanks to the WFC-3 camera on HST, 
provide constraints on the likely sources of hydrogen reion- 
ization. However, these observations only directly reveal 
the sources luminous enough to be individually detected. 
Furthermore, converting the luminosity function measure- 
ments to IGM ionization rates involves large uncertainties, 
in addition to the extrapolation necessary to model sources 
too faint to be detected, including the SED of the star- 
forming galaxies and their escape fraction of ionizing pho- 
tons. These latter uncertainties are encapsulated in the di- 
mensionless factor (^ion/osc used to convert from a luminosity 
at 1,500 A to a rate of production of ionizing photons es- 
caping into the IGM. 

Given these uncertain parameters, we showed in ^that 
many scenarios exist in which star-forming galaxies are the 
dominant ionizing sources and which satisfy both the galaxy 
survey constraints and the Thomson optical depth implied 
by the WMAP-7 data. Such scenarios include ones with 
escape fraction ranging from /esc = 5% to /esc = 50%, 
and limiting UV magnitude ranging from Mum = —16 to 
Afiim = —10, even when these are assumed to be constant. 
Thus, these constraints alone (which have been the focus 
of many analyses; e.g., [Bouwens et al.||2011c| [Bunker et al.| 
20101 are not sufficient to determine the role of faint galax- 



ies in reionizing the Universe, and whether such galaxies are 
even present in significant number. 

In ^we used the Lya forest at redshifts 2 ^ z ^ 6 to 
measure the total instantaneous rate at which ionizing pho- 
tons are injected into the IGM. Although these measure- 
ments cover redshifts past the epoch of reionization, they 
provide significant leverage over galaxy surveys. In particu- 
lar, the total ionizing emissivity implied by the mean trans- 



mission of the Lya forest does not rely on assuming an es- 
cape fraction or a limiting magnitude, two of the main un- 
certainties limiting the predictive power of UV luminosity 
function measurements alone. At z = 4, where the Lya for- 
est data is both abundant and free of large systematic effects 
due to inhomogeneities in the ionizing background, compar- 
ison with the UV luminosity function allowed us to deter- 
mine Cion/esc almost independently of Afum, owing to the 
comparatively shallow faint-end slope of the LF at 2 < 4. 
Since ^ —16 at 2: = 4 (the UV luminosity function 

having already been measured down to this magnitude), 
/,,,(2 = 4) = 2 - 8% (median f^.^z = 4) = 4%), for the 
fiducial spectral model ("ion ~ 1- 

Combining the Lya forest, WMAP-7, and galaxy sur- 
vey data and assuming that galaxies are the main ioniz- 
ing sources requires either: 1) extrapolation of the galaxy 
luminosity function down to very faint UV magnitudes 
Miim ~ —10, corresponding roughly to the UV background 
suppression scale (e.g., Faucher-Giguere et al. 2011; but see 
Dijkstra et al. 2004); 2) an increase of the escape fraction 
by a factor > 10 from 2 = 4 to 2 = 9; or 3) more likely, a 
hybrid solution in which undetected galaxies contribute sig- 
nificantly and the escape fraction increases more modestly. 

The present data do not allow us to select a unique 
viable reionization scenario. Quantitatively, a range of com- 
binations of limiting magnitudes and redshift evolution of 
the parameters affecting the conversion of ~ 1,500 A UV 
luminosity functions to rates of production of ionizing pho- 
tons are allowed and summarized in Figures [7] [8] and |11| 
Redshift evolution in the limiting magnitude alone requires 
extreme assumptions in order to satisfy both the Lya forest 
and WMAP-7 constraints without appealing to very faint 
galaxies. Even so, such scenarios predict that reionization 
ends at 2 < 6, in tension with recent measurements of the ki- 
netic Sunyaev-Zeldovich effect by SPT, which indicate that 



reionization ends earlier than 2 — 5.8 at 95% CL (Zahn et 
ar]|2011 1. On the other hand, significant redshift evolution 
in Cion/csc is more plausible. In fact, Qon can vary by a factor 
~ 4 owing to changes in the age, metallicity, and IMF of the 
stellar populations (corr espo nding to the range (ion = 0.5 — 2 
assumed in this paper; ^2.3| lP]More importantly, /esc can in 
principle increase from /esc ~ 4% at 2 = 4 to /esc ~ 1 at ear- 
lier times. A similar strong redshift evolution of the ionizing 
luminosity-weighted escape fraction was also recently found 
to be required in the "minimal cosmic reionization model" 



of Haardt & Madau (2011 1 



Although there are at present no direct constraints on 
the escape fraction from faint galaxies during the epoch of 
reionization, deep searches for escaping Lyman continuum 
radiation at lower redshifts do show some evidence for red- 



shift evolution ( 


Steidel et al.||2001 


Inoue et al.||2006 


|Shap- 


ley et al.||2006| |Cowie et al.||2009| 


Siana et al.||2010| 


Nestor 


et al. 


2011 


1. Such evolution, in which the escape fraction 



increases with redshift, could owe to increased feedback at 
earlier times when star formation was more vigorous (e.g.. 
Wise & Cen 20091. This picture, in which ionizing pho- 



* Consistency between star formation rate and stellar mass den- 
sity measurements at z ~ 7 — 8 suggests that the IMF of 
reionization-epoch galaxies is not very different from the local 
Universe ( Bouwens et al.|2011b 1 . 
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tons escape galaxies along lines of sight cleared of obscuring 
gas, would be consistent with Lyman continuum observa- 
tions suggesting "on/off" escape, possibly connected to the 
viewing geometry (e.g., Shapley et al. 2006 Nestor et al 



2011 Vanzella et al.|2012 |. Another possibility is that faint 



tion with (-A/um — —16). Thus, experiments capable of mea- 
suring the duration of reionization will have direct implica- 
tions for faint galaxies. In fact, the more extended reioniza- 
tion histories which rely on a heavy contribution from faint 
galaxies are already ruled out by recent constraints from the 
kinetic Sunyaev-Zeldovich effect by SPT. In the next few 
years, expanded data sets from high-resolution microwave 
background experiments in combination with more precise 
measurements of the integrated Thomson optical depth with 
PlancI^wiW improve these constraints further. Refined anal- 
yses should also improve the accuracy with which contam- 
inating point sources are subtracted, and thus solidify the 
results. 

On the theoretical front, there is much room for improv- 
ing our understanding of how to reliably use astrophysical 
sources such as high-redshift Lya emitting galaxies, lumi- 
nous quasars, and 7— ray bursts to measure the neutral frac- 
tion in the IGM. Finally, efforts aimed at detecting 21 cm 
emission from high-redshift intergalactic neutral gas (e.g., 
Dekel| ( |2011| l and [Kuhlen et all ( |2012[ ), metallicity effects IBo wman et al.|2007| [Parsons et al.|2010[ [Bowman fc Rogers] 



galaxies may typically have higher escape fraction than more 
massive galaxies (e.g. Yajima et al.|20rT |, in which case the 
larger relative abundance of faint galaxies at high redshift 
would result in an increase in the population-averaged es- 
cape fraction. The extremely blue UV continuum slopes re- 
cently reported for 2 ~ 7 galaxies (Bouwens et al. 2010b[ ) 
are also suggestive of weak nebular recombination emission, 
which would be consistent with very high escape fractions 
of ionizing photons (but see Dunlop et al. 2011 for a critical 
analysis of the UV continuum slopes). 

Recent models predict that star formation is suppressed 
in low-mass, high-redshift galaxies owing to the metallicity 
dependence of the transition from warm HI to dense molec- 
ular gas. In the fiducial implementations of [Krumholz fc| 



can strongly suppress star formation in reionization-epoch 
galaxies in haloes of mass of Mh ~ 10^ - 10^" M0, corre- 
sponding to Muv ~ —13 to —16 at z ~ 7. Formally, even the 
brighter A/um of these models can satisfy the existing reion- 
ization constraints, but only for strong redshift evolution in 
Cion/eac (Fig. [7| . Siuce such a H2-regulated star formation 
suppression threshold is not far from the current limits of 
HST observations at z ~ 7, deeper integrations have the 
potential to significantly constrain those models. As [Kuhlen] 



et al. ( 2012 1 showed, the exact halo mass below which metal- 



poor dwarfs are suppressed is however sensitive to the de- 
tails of the model implementation. In particular, the relevant 
mass scale depends significantly on the metallicity floor as- 
sumed to model the unresolved effects of metal enrichment 
by early Pop III stars. If the correct halo mass threshold 
is lower by an order of magnitude relative to the fiducial 
([2012 1 and closer to the predictions 
(2011[), i.e. Mh ~ 10® M0 at 2 ~ 7, 



models of Kuhlen et al 
of 



Krumholz & Dekel 



then abundance matching suggests that the turn over in 
the luminosity function would occur instead at Muv ~ — 13 
(see Fig. |3|. For limiting magnitudes in this neighborhood, 
more modest redshift evolution in (^lon/csc can satisfy the 
galaxy survey, WMAP-7, and Lya forest data. Thus, given 
the present implementation uncertainties, H2-regulated star 
formation in high-redshift dwarfs is consistent with galaxies 
reionizing the Universe. 

For each scenario satisfying the observational con- 
straints considered in this work, we evaluated the timing 
and duration of the corresponding reionization history (Fig. 
[9|. Interestingly, the redshift at which the ionized fraction 
reaches 50% is consistent among the different allowed sce- 
narios, Zreion(50%) ~ 10. This redshift is consistent with 
the redshift of instantaneous reionization Zreion = 10.4 ± 1.2 



implied by the WMAP-7 analysis (Komatsu et al. 20111 



On the other hand, the predicted duration of reionization is 
quite sensitive to the fractional contribution of faint galax- 
ies. For instance, our MAX parameterization with A/um — 
— 10 (i.e. with a heavy dwarf contribution) allows scenar- 
ios with 2reion(20%) = 17 aud 2reion(100%) < 4. Ou the 

other hand, models which rely on bright galaxies are pre- 
dicted to reionize the Universe much more sharply, with 
2reion(20%) — Zroion ( 100%) « 3 for our FIT parameteriza- 



|2010 1 are poised to eventually directly map the reionization 
epoch. 
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